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Abstract. We show that the Harnack inequahty for a class of degenerate parabohc quasihnear 
PDE 

dtu — X* Ai{x, t, u, Xu) + B{x, t, u, Xu), 
associated to a system of Lipschitz continuous vector fields X — {Xi, Xm) in in f2 x (0, T) with 
SI C A/ an open subset of a manifold M with control metric d corresponding to X and a measure 
da follows from the basic hypothesis of doubling condition and a weak Poincare inequality. We also 
show that such hypothesis hold for a class of Riemannian metrics collapsing to a sub-Riemannian 
metric limt^o 5e = <7o uniformly in the parameter e > 0. 



1. Introduction 

In two important works, Saloff-Coste |24j and Grigor'yan [TT] established the equivalence be- 
tween a Harnack inequality for weak solutions of a class of linear, smooth coefficients PDE, with 
non-negative symbol, and two key metric-measure properties of the ambient space: (1) a doubling 
inequality for balls in a control metric naturally associated to the operator and (2) a Poincare 
inequality involving a notion of gradient naturally associated to the operator. The Saloff-Coste- 
Grigor'yan results were later extended and applied to a broad range of PDE and variational prob- 
lems (e.g. |18j . and references therein). One of the motivating examples where such theory can be 
applied consists in a class of subelliptic operators associated to a family of Hormander vector fields 
X = {Xi, ...,Xm) and their Carnot-Caratheodory distance •) (see [H], [21] and [25]). 

The purpose of this note is twofold: (a) We show that doubling property and Poincare inequalities 
imply the Harnack inequality for weak solutions of the quasihnear, degenerate parabolic PDE (jl.ip . 
Vice versa, we note that if the Harnack inequality holds for every operator in this class, then by 
virtue of [24], [TT] both the doubling property and the Poincare inequalities hold. Our work 
represents an analogue, in the subelliptic setting, of the work of Aronson and Serrin in [2]. 

(b) We want to show that the constants arising in the doubling and Poincare inequalities are 
stable with respect to a well known and much used, Riemannian approximation scheme, in which 
the Carnot-Caratheodory metric is approximated (in the Gromov-Hausdorff sense) by a family of 
degenerating Riemannian metrics. 

The motivation for our work is to provide the necessary background to study existence and 
regularity for a large class of non-linear degenerate parabolic PDE by approximation with solutions 
of ad-hoc 'regularized' versions of such PDE, in what is occasionally called the vanishing viscosity 
approach. Part (b) shows that the homogenous structure and the constants in the Poicare inequality 
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associated to the geometry underlying the approximating PDE are stable, so that the results in 
part (a) yield uniform estimates on the Holder regularity of the solutions, in the maximum principle 
and in the Harnack inequalities. 

Our proofs rest on a careful extension of ideas and techniques developed in the works [2], [24], 
[llj and [21]. The arguments are technically involved and rest on results spread through a large 
number of papers. Part (a) has been largely developed by the third named author (GR) in his 
PhD dissertation [23]. While completing the final version of the present paper we were made aware 
of recent, pertinent, and very interesting work of Kinnunen, Marola, Miranda and Paronetto |16j 
where an alternative approach to part (a), in a more general setting, is studied. The authors in 
|16j derive Harnack inequalities out of membership in the so-called parabolic De Giorgi classes, in 
the general context of doubling metric measure spaces endowed with a Poincare inequality. It is 
not clear what is the most general setting in which a regularity theory for parabolic PDE would 
make sense and be developed. This problem is intimately connected to finding the most general 
setting for a first order calculus in metric measure spaces [13] . Among other pertinent references we 
mention: the result of Kinnunen and Kuusi [15] where it is proved that doubling and (l,p)-Poincare 
inequality imply the parabolic Harnack inequality for a class of doubly nonlinear equation of p- 
Laplacian type with weights, and the work of Sturm, see for instance [27] where the Saloff-Coste 
and Grigor'yan results are extended to the general setting of local Dirichlet spaces. 

1.1. Analogues of Aronson-Serrin's Harnack inequality and maximum principle. Con- 
sider a distribution of Lipschtiz continuous vector fields X = {Xi, ...,Xm} in a bounded open set 
$7 C W^. We denote by X* their adjoint, and by •) : x (7 — t- M"*" the control distance associ- 
ated to X. In this paper we will assume that d satisfies a doubling conditiorQ ([D]) w.r.t. Lebesgue 
measure and with doubling constant Cd- Correspondingly one finds a number > that acts as 
homogeneous dimension, i.e. volume of metric balls grows at least like a power of A^ of their radius 
(Proposition 12.1( 1)). In this paper we will always assume N > 2. 

We will show that if one further assumes as hypotheses a Poicare' inequality (jPj) with constant 
Cp and the a.e. differentiability of the metric (jLj) with constant Cl, then one has an analogue of 
Aronson and Serrin's results from [2]. Specifically, we establish a maximum principle and a Harnack 
inequality for weak solutions of the quasilinear degenerate parabolic PDE 

(1.1) dtu = X*Ai{x,t,u,Xu) + B{x,t,u,Xu), 
where Ai and B satisfy the structure conditions (j3.2p . 

Theorem 1.1 (Maximum Principle). Let u be a weak solution of (11. ip in the parabolic cylinder 
Q = X (0,T) and assume that there exists M G M such that for every 6 > one has u < M + 6 
in a neighborhood of the parabolic boundary dpQ = Q x t = L) dQ x (0,t). There exists a positive 
constant C depending only on diam{Q), Cd, Cl, Cp, the structure conditions (13. 2p and on Q such 
that for a.e. x £ Q one has 

(1.2) u{x) <M + CK, 
where K = m\ + \\d\)\M\ + {\\f\\ + \\g\\). 

Theorem 1.2 (Harnack Inequality). Let u > be a weak solution of (II. ip in Q and < p < 20R 

with R as in (iDl) . If one has a subcylinder Q-^p = B{x,3p) x (t — 9p'^,t) C Q then there exists a 



See next section for the pertinent definitions 



HARNACK INEQUALITY 



3 



constant C > depending on diam{il), Cd, Cl, Cp, the structure conditions (I3.2p and on p such 
that 

(1.3) maxii < Cmin(n + p^/c), 

Q- Q+ 

where 

Q+ = B{x,p) X {i- p^,t) and Q' = B{x,p) x {t - 8p'^,t- Ip^) 
9 > is defined as in (13. Sp . and we have let k = (||/|| + \ \g\ \ + \ \h\\). 

1.2. Stability of the homogenous structure and of the Poincare inequality in the Rie- 
mannian approximation scheme to a sub-Riemannian metric. We consider examples of 
vector fields satisfying all the assumptions (D), (L) (P). In view of the seminal work in |21j 
and [T^ it is clear that these assumptions are satisfied by smooth distributions of vector fields 



X = {Xi, . . . , Xm) in i?" satisfying the Hormander finite rank condition 

(1.4) rank Lie{Xi, . . . , Xm){x) = n, \/x £ Q,. 
Consequently weak solutions of the equation 

(1.5) dtu = X*Ai{x, t, u, Xu) + B{x, t, u, Xu), 



written in terms of these vector fields satisfy Harnack's inequality and are Holder continuous with 
respect to the Carnot-Caratheodory metric associated to X. 

In many applications (see for instance [22] , [4] , [5] , [8] and [7] ) one is interested in approximating 
the subelliptic PDE with a sequence of elliptic regularizations, so that the corresponding solutions 
converge to viscosity solutions of the original PDE. From the point of view of the metric background 
this approximation corresponds to studying tame approximations of the sub-Riemannian metric 
(and the associated distance function do{-, •)) with a one-parameter family of degenerating Rie- 
mannian metric (associated to distance functions de{-, •)), which converge in the Gromov-Hausdorff 
sense as e — )• to the original one. Such approximations have been studied, for instance, in |19] 
and [20]. We show that for a specific approximation scheme (used in the references above), the 
constants in the doubling and Poincare inequalities can be chosen independently of the parame- 
ter e, thus leading to uniform Holder regularity of weak solutions of the associated approximating 
equations. 

To show a simple, almost naive, example of such approximation as well as a sample application of 
our work, we consider the Heisenberg group with a set of left invariant horizontal vector fields 
Xi,X2 such that Xi,X2,X2; = [Xi,X2] generates the whole Lie algebra (see [26] for a detailed 
description) . For e > we consider a family of left-invariant Riemannian metrics in defined 
so that the frame Xi, X2,eX3 is orthonormal. A well known technique to study solutions of the 
degenerate parabolic PDE (to mention one example out of many) 

(1.6) dtu = Xfu + X^u 

is to consider a family of solutions u*^ of the strongly parabolic PDE 

(1.7) dtu^ = Xfu^ + X^u" + e^Xfu' 

with fixed (e— independent) data on the boundary of given a bounded parabolic cylinder in Hlix]R+. 
The PDE (11. 7p is the heat equation for the Riemannian metric ge, while (11. 6p is the subelliptic heat 
equation for the sub-Riemannian metric go. Parabolic theory yields existence and uniqueness of 
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smooth solutions u*^ to ()1.7p . The problem is that such smoothness may degenerate as e — )• 0. In 
this paper we show (see Proposition 14. 4p that the doubling property and the Poincare inequality 
corresponding to the Riemannian metrics are stable as e — )• and as a consequence of the stable 
Harnack inequality Theorem 1 1 . 2 1 the set {u'^}e>o is a normal family, with uniform limit uq. Applying 
the same argument to right invariant derivatives of the solutions (see for instance [3] ) we can easily 
prove that uq is smooth and solves (11. 6p . Obviously, the linear structure of (11. 6p provides far more 
effective ways of proving existence of smooth solutions, however such methods cannot be applied 
for non-linear PDE, while the techniques in the present paper are purposedly designed to be used 
in the non-linear setting. 

The core of argument in the stability proof consists in a careful extension of the arguments 
in [21] to include the additional parameter e, and show e— independent bounds on the Jacobian 
of an exponential mapping related to commutators of the Hormander vector fields. Once such 
bounds are established, the doubling condition follows immediately and the Poincare' inequality 
is a consequence of the work in [T7] . To better appreciate the stability results, we note that even 
in the simple example described above the geometry of the space changes radically as e — ?■ 0. In 
this approximation the curvatures tensors blow up, and the shape and volume of the metric balls 
changes from Riemannian to sub- Riemannian. Yet, the homogeneous structure, as well as the 
constants associated to the Poincare inequality remain the same. 

Acknowledgements. The authors would like to thank Juha Kinnunen, Niko Marola and Michele 
Miranda Jr. for sharing with us helpful comments and pertinent references. 

2. Preliminaries 

Let Xi, Xm be Lipschtiz continuous vector fields in a open set C M". We denote by 
d{-, •) : $7 X — )• M"*" the control distance associated to X and for x S and r > 0, by B{x, r) = 
{y £ M"'\d{x,y) < r} the corresponding metric balls and by \B{x,r)\ their Lebesgue measure. The 
average of a function u over a ball B = B(x,r) is denoted by ub- For 1 < p < oo we define the 
Sobolev spaces W^'P{n) = {u £ LP(0)| Xu = {Xiu, ...,X^u) £ LP{n)} and Wq'^ to be the closure 
of C°°{Q) in the norm = \ \u\\p + 

In our results we will rely on the following hypothesis on the background geometry induced by 
X: For every compact K C O there exists constants Cd = Cd{X, K),Cl = Cl{X, K),Cp = 
Cp{X, K) > 0, and R = R{X, K) > such that for every x £ K and < r < R one has 

(D) \B{x,r)\>CD\Bix,2r)\. 

(L) d{-,x) is differentiable a.e. in and | |Xd(-, x)| |ioo(^-) < Cl for every x £ K. 

(P) f \u-UB\^dx<Cp f \Xu\'^dx. 

JB(x,r) JB{x,2t) 

In the following we say that a constant depends on Co-,Cl and Cp if it is uniformly bounded 
when both the parameter R and such constants are so. As a consequence of these hypotheses one 
has that the metric measure space (X, d, dx) is an homogenous space with a Poincare' inequality, 
and as such it enjoys several useful properties listed below. 
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Proposition 2.1. If hypotheses (jD]),(jL]), ([P]) hold then for every compact subset K C Cl, x £ K 
and < r < R one has: 

• (i) (Lower hound on volume of halls) There exists a constant N = N{Cd) > 0, called 
homogeneous dimension of K w.r.t. {X,d,dx), such that \B{x,r)\ > Cgt^ ■ 

• (a) There exists a test function cj) G Wq '°°(i?(x, 2r)) and a constant C = C{Cl) > 0, such 
that (p = 1 in B{x,r) and \X(j)\ < C/r. 

• (Hi) (Weighted Poincare inequality, [ISl Theoreme 8.1]j Let (f) : — t- [0,1] be a non 
increasing function with compact support in a finite interval [0, R) and such that (1) t/0 G 
C^([0, R)); (2) For < r < R one has 4>{r + i min(i? - r, §)) > a(j){r). Then there exists 
C = C{Cd,Cp, Cl,N,R)>0 such that 

(WP) j \uiy) - u^l x))dy < CR-'^/^i j (t^{d{y, x))dyf'^ j \XuWy)(l,{d{y, x))dy, 

for every function u W^''^{Q,) and where we have let K = infy^K \B{y, R)\{2R)~^°^'2'-^'^ >0 

• (iv) (Sobolev Embedding [I8l[l2]j There exists C = C{Cd,Cp, N) > such that for every 
l<p<N, 

N~p 1 

(S) { / \u\^dx) <Cr( [ iXul^dx] 

\JB(x,r) J \JB{x,r) J 

for every function u G Wq'^{Q,). 

We also notice that in view of Cheeger's Rademacher's theorem [6], hypothesis (jLj holds as a 
consequence of ([D]) and (jP]) . 

Next, we consider a parabolic cylinder Q = Q x {0,T) and address the parabolic BMO spaces 
and the John-Nirenberg lemma for the homogenous space {Q,dp,dx), where dp{{x , t) , (x , t)) = 
max{d{x, x), -s/ft^-tl}. 

For (x, t) £ Q and r > sufficiently small we consider a cylinder 

Qrix, t) = B{x, r)x{t- 8r2, t) C Bp{{x, t), lOOr) := {{y, s)\dp{{x, t), {y, s)) < lOOr} C Q 
and define the upper and lower cylinders 

(2.1) Q;{x,t) = B{x,r) X {t - 8r^ ,t - 7r^), {x,t) = B,{x,r/2) x {t - ^r^ ,to - 7r^), 

Qt{x, t) = B{x, r)x {t- r'^,t) and D:^{x, t) = B{x, r/2) x {t-r'^,t- ^r^). 

We define the time lag function T{{x,t),r) = {{x,t — 8r^),r) and observe that for all {x,t) G Q 
one has Q~{x,t) = T{Q^{x,t),r)^ thus satisfying [H Hypothesis (1.8)]. In view of the doubling 
hypothesis (jD]) one can apply [U Theorem 1.7] and conclude the following version of the John- 
Nirenberg lemma 

Proposition 2.2. If there exists A > andv G L^{Q), such that for all cylinders Q2r,Qr ^Qt C Q 
as defined above there is a number C = C{Q) such that 



(2.2) — — / / ^{v - C{Q))+dxdt < A and — — / / ^/(C(^^)^^ dxdt < A 

\Q \ Jq+ j \Q I Jq- 
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then there exist C = C{Q,Cd,^),S = S{^1,Cd, ^) > such that if we set f = e~^ then for all 
D^,D~ C Qsr C Q as defined above 

3. Maximum principle and Harnack inequality 

Let C M" be a bounded open set and let Q = 0, x (0,T). For a function u : Q — )■ M, and 
1 < p,q we define the norms 

(3.1) \\u\\'^p,q = [ if \u\Pdx)pdt, 

Jo Jq 

and the corresponding Lebesgue spaces L^''^{Q) = L'^([0, T], LP(ri)). To simplify notations we will 
omit the subscripts p, q when their values are clear in view of the context. One has a useful refor- 
mulation of the Sobolev embedding theorem in terms of L^'"? spaces and a related interpolation 
inequality: 

Lemma 3.1. (i) Let Xu G L?'''^{Q) and assume that for all < t < T , u{-,t) has compact support 
in Q X {t}. There exists K = K{diam il.,CD,Cp, N) > such that \\u\\ 2n n < -fC||Xn||22- Here 

N is the homogenous dimension from Proposition \2. l]( i). 

(a) If in addition we assume u E L'^'°°{Q) then u G L'^p ^"^Q (Q) where {p',q') have as Holder 
conjugates {p,q) such that N /2p + \/q < 1. Moreover there exists K > depending only on Cd,Cp 
and diam Q such that 

\\n\\%',2,' < KT%\\w\\l^ + \\Xw\\l^), 

where 9 = 1 - q'^ - N/2p. 

We will consider weak solutions of (jl.ip . where Ai and B satisfy the following structure condi- 
tions: There exist constants a,a > and functions b, c, e, f,h & LP''^{Q) with p > 2, and q given by 
^ + I < ^ and functions d,g €z L^'^iQ) with 1 < a and /3 given by ^ + ^ < 1 such that for a.e. 
(x, t) £ Q and ^ G one has 



^A{x,t,u,0^^ >a|ep-6V-/2 



i=l 

(3.2) \A{x,t,u,^)\ < a\^\ + e\u\ + h, 

\B{x,t,u,0\ <c\C\+d\u\+g. 
In view of the conditions on p, q, a, /? there exists 9 > such that 

2 n 1 1-9 

p > ^ and - — I — < 



1-9 2p q- 2 

(3.3) and^ + l<l-.. 

In this note when we say that a constant depends on the structure conditions (|3.2p . if it depends 
only on 

a,a,||6||,||c||,||d||,||e||,||/||,||5||,||/i||,A^,e, 



The |[ ■ II norms are in the appropriate L*"'' or L"''^ classes 
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and is uniformly bounded if these quantities are so. Here N denotes the homogenous dimension in 
Proposition I2.1( i) . 

A function u G L'^'°°{Q) with Xu G L^'^((5) is a weak solution of (ll.ip in Q if 



(3.4) 



-u(j)t + Xi(j)Ai{x, t, u{x, t),Xu{x, t))dxdt 



4>{x, t)B{x, t, u{x, t),Xu{x, t))dxdt, 



for every (p G Cq°{Q). To deal with the lack of differentiability of weak solutions along the time 
variable we will make use of Steklov average: For (x, t) £ Q and h > sufficiently small set 
Uh{x, t) = Jq u{x, t + s)ds. Changing variables t — >• t + s and integrating in s in (13. 4p leads to 



(3.5) 



dxdt 



dt{uh)4> + Xicph ^ I Aiyx,t — s,u{x,t — s),Xu{x,t — s)jds 

J J (f>{x,t)h~^ J B^x,t — s,u{x,t — s),Xu{x,t — s)^dxdt, 



for every cj) G C^{Q). Clearly the identity extends to the larger class of test functions with weak 
derivatives in L^'^ and with vanishing trace. 

Lemma 3.2. Let u be a weak solution of (jl.ip in Q = Q x (0, T) and k a positive constant. For 
every r] G C°°{Q) vanishing in a neighborhood of the parabolic boundary dpQ, one has 

• If P u = max(0, u) + K and < t < T , then 



(3.6) 



1 



P + 1 



7] 



t=T 



dx + - 



Jn 



rj^uf^^^lXul^dxdt 



< 



Jn 



Fu^^^dxdt + 



/3 + iyo 



ri\dtrj\u^^^ dxdt., 



where T = Frf + 2Gri\Xr]\ + H\Xri\^ , and 



F = /3( 6^ + ^ ) + ( d+ii ) +-, G = e+-, andH = —. 



c 
a 



h 



4a^ 
a 



If P ^ 1) u = max(0, u) + K and < ti < t2 < T then 



(3.7) 



u^^^ - (/3 + l)k'^u + /3A:^+i 



''dx+^ r [ r^'uf'+^lXul^dxdt 
ti 2 Jti Jn 



< 



ti Jn 



Tu^+^dxdt + 



t2 



ri\dtri\u^^^dxdt, 



with T defined as above. 

If u > and for arbitrary e>0 we let u = u + k + e, and < ti < t2 < T then 



(3.8) sign{f3) 



ifn{u)\\idx+''4- 



t2 



j]'^u^-^\Xu\^dxdt 



< 



t2 



u Jn 



Tiul^^^dxdt + 2 



t2 



u Jn 



ri\dtri\\7i{u)\dxdt, 
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where 



(3.9) 
(3.10) 



niu] 



^^7/3+1 



^//5 = -l, 



I log u 

Fir?2 + 2Gi77|Xr/|+//i|Xr/|2, 



with 



Fi 



9 



Gi = e + — , and Hi 

K 



4a' 



Proof. The argument follows closely [21 Section 2] with Steklov averages in place of the elegant 
convolution argument in Aronson and Serrin's original paper, as the latter does not go through as 
easily in the case of variable coefficients vector fields. Let 

— for — oo<u<l — K, 

l^~^u — for I — K < u < oo, 



Giu) 



where I > k. Denote by s — )• H{s) a function such that H'{s) = G{s) for all s G M. Choose 
r] G C°°(Q) vanishing in a neighborhood of dpQ and set 4>{x,t) = 'rf'{x,t)G{uh{x,t))x{o,T){t) in 
), to obtain 



(3.11) 



if'H{uh)dx 



t=T 



Jn 

h 



2rjdt'q'H{uh)dxdt+ 



I I Xi(f)h ^ / Ail x,t — s,u{x,t — s),Xu{x,t — s) jds 
Jo Jn Jo \ J 



dxdt 



j j (f){x,t)h ^ j B^x,t — s,u{x,t — s),Xu{x,t — s)^dxdt. 



At this point we can let /i — )• and make use of the structure conditions to conclude (13. 6p in the 
same fashion as in [2]. The estimates (13. 7p and (13. Sp are proved in the same way. □ 

3.1. Proof of the maximum principle. The proof of Theorem ll.21 is accomplished in two steps: 

Lemma 3.3. If M < then there exists a positive constant Gi depending only on Gd,Gl,Gp and 
the structure conditions (j3.2p such that the function u = max(0, u) satisfies 

(3.12) u{x,t) <Ci(||u||2,oo + ||^u||2,2 + k), 

for a.e. {x,t) G Q and with k = ||/|| + ||(7||. 



Proof. As in proof of Lemma 13.21 we set u = u + k and note that since Giu) vanishes in a neigh- 
borhood of the parabolic boundary of Q then (|3.6p holds with no need of the cut-off function 
V, 



(3.13) 



1 



/3 + 1 



dx + 



t=T 







f+^\Xu\'^dxdt 



< 



Jn 



Fu'^^^dxdt 



'The same dependance holds for the constants C2, C3, ... used in the proof 
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where F = (^(^b'^ + + (^d + + Set r = (/3 + l)/2 and v = vJ' . Letting t ^ t \n (l3J3]l 

one obtains a| |2,2 — i. Let |||it;||| = supp^ ||w^||p,g where the sup is taken over all the 

pairs {p,q) satisfying n/2p + 1/g < 1 — 9 and p > 1/(1 — 0). In view of the structure conditions 
(j3.2p one has that there exists C2 > such that /3~^||F?;^||i^i < C2|||t^||p, and consequently 

(3.14) a||Xi;||^^2 < <^2r^||blll^. 

Young inequality and p.lSp yield ||f||2oo ^ C'3r^|||t;||p for some C3 > 0. The latter, Lemma [XT] 
and ()3.14p imply 

(3.15) \\\v\\r <m+m\v\\loo + \\xv\\y<cy\\\v\\\^ 

with £7 = 1 + 26 /N. A standard iteration argumenlQ in the exponent of v leads to (I3.12p . □ 

The second step concludes the proof of the maximum principle. 

Lemma 3.4. If M < then there exists a positive constant C depending only on Cd,Cl,Cp and 
the structure conditions (|3.2p such that the function u = max(0, u) satisfies 

(3.16) \\u\\l^ < Ck^ and \\Xu\\\2 < Ck^, 
with K = 11/11 + ll^ll. 

Proof. Let /3 = 1 in (I3.7p . with the cut-off function rj omitted, and apply Holder inequality and the 
interpolation in Lemma 13.1( 11) to obtain 



dx -\ — / / \Xu\ dxdt 
ti 2 Jti Jn 



ft2 

< I I Fu^dxdt 



ti Jn 



< 3iK + l)C{t, - t^YiMl^ + \\Xu\\l, + k'), 

with 6 as in (j3.3p and for a constant C as in the statement. Next, we let fi = [min(a, 1){K + 
l)C/8]^/^, choose ti < t < t2 + fJ. and set X{t) = j^v?{x,t)dx. Estimate (|3.17p yields 

(3.18) X{t) + ^ f [ \Xu\'dxdt < (||72||2,oo + K^) + X{t,). 

^ Jti Jn ^ 

Iterating this inequality from one time interval to the next yields Xii) < 2^~^^^ ^ which for t = T 
gives the first inequality in (j3.16p . The second inequality follows immediately from the latter and 
from (f3l7D . □ 



^See [2] page 95] for details 



10 LUCA CAPOGNA, GIOVANNA CITTL AND GARRETT REA 

3.2. Proof of the Harnack inequality. We let Q,., Q^,Df^ denote the cylinders defined in ()2.ip . 

Lemma 3.5. Let u > be a weak solution of (|l.ip in Q = Q x (0, T), 9 be as in (j3.3p . the 
homogenous dimension (Proposition \2.lV i)) and < p < 20 R with R as in ([D]1 . There exist a 
constant C > depending on diam{Vt) ,C o , C p , C l and on the structure constants (j3.2p such that 
for every Q2p{xo,to) C Q and every exponent of the form tq = (1 + 2-^)~''(2 + 2-^)^-'^, h = 1,2, ... 
one has 

(3.19) esssupjy-\u^°\ < C sup\\if°\\^pi,qi ^Q--^ = ClllvJ' \\\q- , 

§ p',q' ^ § 

and 

(3.20) ( essinfjj+\u^'^\] < Csup ||'u"'""||^p/,,/(Q+) = C|||n~'^°|||Q+ , 
\ f / p',g' f f 

where the supremum on the RHS is taken over all pairs {p',q') whose Holder conjugate satisfy 
n/2p + 1/q < 1 — 9 and p > 1/(1 — 9), and where we have let u = u + k + e with e > arbitrary 
and k as in the statement of Theorem M.Si 



Proof. Since translations in the time variable do not affect the structure constants in (|3.2p we will 
prove the two estimates (j3.19p and (|3.20p separately in cylinders 

f 1 — s 1 + s 
S{s) = B{xo, sp) X -^P^ 



^ 6 ' 6 

for each 1/3 < s < 1/2 Let /3 / -1. Substitute in ([M]) v = u'' with r = (/3 + l)/2 to obtain 

ft2 



(3.21) sign(/3) 



/ r]^v'^\l^dx + ^ /' / rj^lXvl^dxdt 
Jn 2r^ Jt^ Jq 



< r [ {Ti + 2\P + 1\-Wdtr]\ydxdt. 
Ju Jn 



/tl Jn 

To prove (pM) we let /? > -1. For 1/2 </'</< 1/3 choose rj G Lipo{S{l)) with 77 = 1 on S{1') 
and \Xri\ < (i^ii-^p and \dtri\ < jp^r^- Substituting 77 in p.2ip and using Holder inequality we 
obtain the following estimate for the RHS, 

(3.22) \\{Ti+2\P + l\-Wdtr]\y\\i,i 

< max(|/3|, 1/3 + \I3\)C{1 - /')" 

where here and for the rest of the proof we indicate by C constants as in the statement of the 
lemma. In view of p.22p we obtain for all /3 > — 1, /3 7^ 



C(l + |/3|-2 + |/3|-i)(l+)fl + (^)^ 



2-r 

(3.23) ||??^'^||2,oo + a|k^^^ll2,2 < pi{l-VY ^ ^lll^lllsw 

By means of the Sobolev inequality (Lemma 13. ip the latter yields the basic iteration formula 



(3.24) |||(^<||p/^< 



C(l + |/3|-2 + |/3|-i) l + (^)2 



2 
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where the gain in integrabihty is a = 1 + 2-^, and N is the homogenous dimension. Set the iteration 
step to be = roa"^, for m = 1,2, .... To avoid the exponent /3 = 0, corresponding to r = 1/2 
it is convenient to set r^^ = a^{l + cr) for any h = 1,2,.... To carry out the iteration we let 

= I + — and = i + — and proceed as in Moser's original paper arriving eventually 
at the estimate 

esssup5(i)|u'^«| < C\\\u''° Wisely 
The proof of (|3.2U|) is very similar and is omitted. □ 

Next, we turn our attention to the 'bridge' step in the Moser iteration scheme, where the sup 
and the inf estimates established above are linked by means of the John-Nirenberg lemma. 

Lemma 3.6. Let u > be a weak solution of p.l|) in Q = Q x (0, T) and < r < 20i? with R 
as in ([D]) . There exists B > depending on diam{Q),CD,Cp,CL and on the structure constants 
()3.2I) such that for every Q2rixo,to),Qrixo,to),Q~{xo,tQ) C Q one has 



(3.25) TTTfTTT^ / / \ ( '^og u{y, s) - log u{x,t)] dyds dxdt < B, 

\Qr\ \Qr I jQ+rJQ- V V / 

where we have let u = u + k + e with e > arbitrary and k as in the statement of Theorem 
Proof. Set /3 = — 1 and v = log u in ()3.8p to obtain 



(3.26) - / if{v{-,t2)-v{-,ti))dx + - / / 'q^\XvYdxdt 

'■t2 r r*2 



n 2 jQ 



t2 

2| 



rt2 r rt2 r 

< / Fidxdt + 2 / / 7]\dtri\\v\dxdt 



h 

where Fi is as in (|3.8p . Next we choose the cut-off function so as to match the requirements of 
the weighted Poincare inequality (jWP]): r]{x,t) = ^{x)w{t) with ^(x) = \/(j){d{-, xq)) and where (j) 
satisfies the hypothesis in Proposition I2.ir iii). In this way ^ E LipQ{B(xo,r)) and can be chosen 
identically equal to one in B{xQ,r /2). We also let w = 1 if t > ti and w = for t < ti/2. Setting 

_ J^^Mx,t)dx 

in (j3.26p and applying (jWPp one obtains 

n ck"^/'"- /"^^ r 

(3.27) V{t2)-V{t^) + - \v-V{t)\Ydxdt 

^ / \ N Jti Jn 

< pi, / / Tidxdt, 
in^^dxJt, Jn 

wherein = inf^ef^ |5(x, r)|r-'°S2 > q. Since \B{xo,r)\-^/^ > C{diam9) > the estimate ([3:271) 
yields 

(3.28) ^ + / {v- Vfdx < — ^ / T,dx, 

dt \B{x,r)\ JB(xo,r) \B{x,r)\ JB{x,r) 
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for every < t < T and for constants Ci,C2 > depending on the structure conditions and 
on diam{Q),CD,CL and Cp. The conclusion now follows from [21 Lemma 7] as in the Euclidean 
setting. □ 

The previous lemma tells us that logn is in the parabolic BMO space. Applying Proposition 12.21 
we conclude 

Corollary 3.7. Let u > be a weak solution of (jl.ip in Q = x (0,r) and < r < 20i? with R 
as in (jD]). There exist C, (5 > depending on diam{^) ,C d ■, C p , C l and on the structure constants 
()3.2p such that for every Q2r{xo,tQ),Q:^{xo,tQ),Q~{xo,tQ) C Q one has 



(3.2.) (_^/^ Je-...,) < c. 

where we have let u = u + k + e with e > arbitrary and k as in the statement of Theorem \1.SX 

In order to complete the proof of the Harnack inequality one needs to link ()3.29p with the RHS 
of inequalities (I3.19P and (I3.20p . This can be accomplished following the same argument at the end 
of page 106 in [2]. 

4. Stability of the homogenous structure and of the Poincare inequality in the 

RiEMANNIAN APPROXIMATION TO A CaRNOT-CaRATHEODORY SPACE 

Let X = (Xi, Xm) be smooth vector fields in M" such that, together with all their commutators 
up to step r, they generate the whole M". Following |2H page 104] we define the collections of 
commutators of same degree 

X« = {Xi,...,X^}, ={[Xi,X2],...,[X^_i,X„,]}, etc.... 

Indicate by Yi, ...,Yp an enumeration of the components of X^^\X^'^\ ...,X^'^^ such that Yi = Xi 
for every i < m. If Yfc G X^'''^ we set the degree of Yk to be d{Yk) = d{k) = i. 

For each 1 > e > e > consider rescaled vector fields, and a suitable subfamily of their commu- 
tators 



y^ ifi<m, I Xt ifi<p, 

I \ e'i(0-iy. iim + l<i<p' * | ^i-p+m if p + 1 < i < 2p 



m 



We will also extend the degree function, setting d^ii) = 1 for all i < p, and d^{i) = d{i — p + m) if 
i>p+l. In order to simplify notations we will denote X = X^ , Y = y*^, d^ = d and use the same 
notation for both families of vector fields (dependent or independent of e). 

Note that for every e the sets {Y^} extends the family (X|) to a new family of vector fields 
satisfying assumption (I) on page 107 [21]: There exist smooth functions c^-^, depending on e, such 
that 

[Y-,Y^]= YI 4fc^ 

d4l)<d4j)+d^{k) 

and 

{Yj}'^^"^ span M."' at every point . 

Remark 4.1. Note that the coefficients c^^ will be unbounded as e — )• 0. In principle this is a 
problem for us, as the doubling constant in [21] depends indirectly from the norm of these 
functions. 
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Next we consider the Carnot-Caratheodory metric de{-, ■) associated to the family of vector fields 
{Xi, ...,Xp). Note that for 1 > e > e > 0, tie corresponds to the distance function of a Riemannian 
metric g^, which degenerates to a sub- Riemannian metric as e — )• 0. 

Following (21] page 110], for every n— tuple / = (zi, ...,in)i we define the coefficient 

XKx) = detiY^{x),...,Yl{x)), 
for e > 0. For a fixed constant < C2,e < 1, choose = {i^i, ...,ien) such that 

(4.1) |A^^(x)|r'^'(^^) > C2,,maxj|A}(x)|/^(-^). 

and denote the family of remaining indices, so that {Yi^j '■ iej G -^t} U {i^i^^. : iek G Je} is the 
complete list V". We will refer to Iq as the choice corresponding to the |21j setting. The core of 
the result of Nagel, Stein and Wainger |2H Theorem 1], is to prove that if v and x are fixed, and 

Qe{r) = {u G : \uj\ < r'^^fe)} 

is a weighted cube in M", then |Aj^(x)| provides an estimates of the Jacobian of the exponential 
mapping — )• ^e,v,x{u) defined for u £ Q{r) as 

(4.2) $.,.,.(«) = exp( J2 ^jK, + E ^fc^4.)(^)- 

More precisely, for e > and fixed |21l Theorem 7] states 

Theorem 4.2. For every e > 0, and K CC M" there exist > and constants < Ci^^, C2,e < 1 
such that for every x £ K and < r < R^, if is such that (j4.ip holds, then 

i) if \vk\ < C2er'^^*"=\ ^e,v,x is one to one on the box Q^iCi^^r) 

ii) if \vk\ < C2er'^^*'=''^ the Jacobian matrix of^e,v,x satisfies on the cube Qe{Ci^ef) 

\\Xlix)\<\J<^>,^,^,\<A\\l(x)\ 

iii) 

^e,v,x{Qe{Cl,er)) C B,{x,r) C ^e,vAQe{Cl,er/C2,e)) 

A direct consequence of this fact is that the measure of the Ball centered in x can be estimated 
by the measure of the cube and the Jacobian determinant of ^e,v,x- Varying the central point x, 
[21] obtain: 

Theorem 4.3. f |21l Theorem 1]) For every e > 0, and K CC M" there exists a constant R^ > 
and constants Cse, C^e > such that x £ K and < r < R^ 

(4.3) Cs. j;|AK^)|r'^(') < mx,r)\ < |AKx)|r'^(^), 

To prove that that these inequalities hold uniformly in e, it is enough to study the constants Cie 
C2e and show that they do not vanish as e — )• 0. Without loss of generality we can assume that both 
constants are non-decreasing in e, otherwise we consider a new pair of constants Cj^^ = inf^g[^ Cj^^, 
fori = 1,2. 

Proposition 4.4. For every e G [0,e], the constants i?e)C'i,e (^nd C2,e in- Theorem \4-^ may he 
chosen to he independent of e, depending only on the C^~^^ norm of the vector fields, e, and on K . 
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Proof. The proof is split in two cases: First we study tlie range e < r < i?o which roughly 
corresponds to the balls of radius r having a sub-Riemannian shape. In this range we show that 
one can select the constants to be approximately Cifi. The second case consists in the analysis 
of the range r < e < e. In this regime the balls are roughly of Euclidean shape and we show that 
the constants Ci^^ can be approximately chosen to be Ci^i. 

Let us fix e > 0, i? = i2o ^-nd r < Rq. We can start by describing the family defined in (j4.ip . 
which maximize X}{x). We first note that for every e > and for every i, m + 1 < i < p we have 

(4.4) y/r^^^^*) = e'^^^~\Yi, y/+p_^/'(^+P-") = r^{i)Y,. 

Case 1: For every e < r < Rq the indices 1^ defined by the maximality condition (14.10 coincide 
with indices of the family Iq and do not depend on e. On the other hand 

J^ = jQU{k :m+l <k<p} 

so that 

{Y^-.ke JJ = : io,k G Jo}U{e'^(*o-^)-iy,o , : io,k G Io,io,k > m+l}U{e'^(*°-^)-^yi„ , : io,k G Jo}- 

In correspondence with this decomposition of the set of indices we define a splitting in the u— variables 
in ([42]) as 

V = {v, V, v). 

Consequently for every e < r the function ^e,v,x{u) reduces to 

(4.5) <I>,,,,,(^z) = exp( ^jK, + E ^fc^4)(^) = ^^^'( E ""K, + E ^^^^4)< 



Let us define mappings 
id 



In view of ()4.5p we can write 

(4.6) ^e,v,xiu) = %,F2Av)APh^Au))- 

Note that for any e > and for a fixed v, the mapping u — )• Fi^^^i,{u) is invertible and volume 
preserving in all M". Moreover J$e^„^a;(u) = J^Q p^ ,j.(i^i^e^t,(u)). In view of (|4.6p and of Theorem 



21 as a function of u, the mapping ^e,v,x{u) is defined, invertible, and satisfies the Jacobian 
estimates in Theorem 14.21 (ii) 

\\Xl{x)\ < \J%,F,Av)APlM^))\ = \J^e,vA^)\ < 4|A?„(X)| 

for all u such that Fi^^^^{u) £ Qo{Cifir) and for v such that 
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|ni| < Ci,or'^(*°i) • • • \um\ < Ci,y^'^^\ \um+i + Wie''^''""+^^"'l < Ci,o/(*0'"+^\ 
when /c = 1, 2p — m. 

The completion of the proof of Case 1 rests on the following two claims: 

Claim 1 let e < r < Rq. There exists Cq > 0, independent of e, such that for all v satisfying 
\vk\ < Cgr^*'-'') one has \F^^^{v)\ = \vk + Vke'^^''"'^~^ < C2,or'^^'°'<^l 
Proof of the claim: If we choose Cg < minjCi^Oj ^^2,0} and 



it follows that 



l-Vfcl, bfel < niin{Ci,o,C2,o} — ^ — , \uj\ < Ci^- ^ 

\Vk\ < C2,0 \Vkl \vk\ < Cl,0|, < Cifl ^ . 



So that 



\vk\<C2,o^^, e''^'°'^-'\vk\, e'^('«^)->fc| < Ci,o^^, \uj\ < Ci,o^^ 



completing the proof of the claim. 

Claim 2 Let e < r < Rq and v fixed such that \vk\ < Cqt'^^'^^^^ for k = 1, 2p — m. One has 
that 

Q,(C5-V) c Fi-;,(Qo(Ci,or)) C Q,{C^t) 

for some constant C5 > independent of e > 0. 

Proof of the claim: Choose C5 sufficiently large so that 2 max{ C^^jCe} < Ci^ and observe 
that if n G Q^{C^^r) then for k = l,...,m we have \uk\ < Ci^r'^^^'^''^^ = Ci^or'^^'''0-'''> while for k = 
m + 1, ...,n we have \Fl^^^{u)\ = litfc + 'Sfce'^(*ofe)-i| < max{C^^ ,CQ}r'^^^ok) (^i^^d{iok)-i^ < Cifir'^^'o''K 
This proves the first inclusion in the claim. To establish the second inclusion we choose C5 large 
enough so that 2(Ci^o + C'2,e) < C's and observe that if Fi^^^i,{u) £ Qo{Ci^or) then for k = m + 1, ...,n 
one has \uk\ < |«fc + i;fce'^(*f"^)"^| + |i;fc|e'^(*o'=)"^ < 2{Ci^o + C2,e)r'^^'°''^ < C^r'^'<'°^\ The corresponding 
estimate for the range k = 1, ...,m is immediate. 

In view of Claims 1 and 2, and of Theorem 14.21 It follows that for e < r and these choice of 
constants (independent of ejlthe function is invertible on Qo{Cifir) and i), ii) and iii) are 

satisfied. 

Case 2: If r is such that r < e < e we can apply and deduce that, if 

0^ )i0,l&Io ~ {-^0,1' ■ ■ ■ !^0,n} 

then the maximality condition (|4.ip selects the following family of vector fields: 

('V'^V — Sf:diio,l)-lY. ... ,rf(io,n)-lA^. 1 

The complementary family becomes 

: Uk G Je} = {e'^^'"'''^~%,,, : io,k G Jo} U {Yi^ , : io,k € Jo} U : io.fc e /o,io,fc > m + 1} 

If we denote A^, and Ce these three sets, and split the f— variable as u = ({), then it is 
clear that 

fd{io,k)-l 

y G A, iff . Y G A^, 



^Ro in place of R^, C5 in place of Ci,e and Cg in place of C2,e 
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and in this case the values of and di are the same on the corresponding indices. Analogously 

and the degrees are the same. 

Consequently in this case, for every e > r the function ^e,v,x{u) reduces to 

#,,,,,(n) = exp( J] u,Y^^+ VkY{l^){x) = expi^ ^ ^^.e4io,.)-iy.o^ + ^ VkY^^){x) = 

/ grf{*0,fc)-l _ gd(io,fc)-l _ s 

M E ^^-^^(^^^4 + E + E ^^k + E ^^^^^(^^4)(-) 

This function is defined and invertible for 

Recall that with the present choice of r < e < e, we have Ci^e?"**-*"^^ = Ci^^r'^^^^^^^ = Ci^ir'^'^'^^^\ If 
we set 

and argue similarly to Case 1, then the function ^e,v,x will satisfy conditions i), ii), and iii) on 
Q{Ci^ir) and hence on (^(Ci^^r), with constants independent of e. □ 



We now turn our attention to the Poincare inequality and prove that it holds with constant 
independent of e. Our argument rests on a rather direct proof from [17j which in some respects 
simplifies the method used by Jerison in [13] • Using some Jacobian estimates from |10j or [9] 
we will establish that the assumptions required in the key result [171 Theorem 2.1] are satisfied 
independently from e > 0. We start by recalling 

Theorem 4.5. [I7i Theorem 2.1] Assume that the doubling condition (D) is satisfied and there 
exist a sphere B^{xq, r), a cube C M" and a map E : B^{xQ,r) x — )• satisfying the following 
conditions: 

i) B^{xo,2r) C E{x,Q^) for every x G B^{xo,r) 

ii) the function u i— )■ E{x, u) is one to one on the box Qe as a function of the variable u and 
there exists a constant ai > such that 

— \JE{x,Q)\ < \ JE{x,u)\ < ail J^(x,0)| for every u e Qe 

Also assume that there exists a positive constant a2, and a function 7 : B^{xQ,r) xQ^x [0, 02?"] — t- 
satisfying the following conditions 

iii) For every {x,u) £ Be{xo,r) x the function 1 1— )• j{x,u,t) is a subunit path connecting x 
and E{x, u) 

iv) For every {h,t) £ B^{xQ,r) x the function x 1— )• ^{x,u,t) is a one-to-one map and there 
exists a constant 03 > such that 



inf 

Be{xo,r)xQ^ 



det-— 
ox 



> as 
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Then there exists a constant Cp depending only on the constants ai,a2, as and the doubling constant 
Cd such that (P) is satisfied. 

Proposition 4.6. The vector fields {Xl)i=i...p satisfy condition (P) with a constant independent 
ofe. 

Proof. We use Proposition 14.41 to show that the assumptions of Theorem 14.51 are satisfied unformly 
in e. Fix a value e > and a set K = Bf^(xo,r). Choose the constants Cj as in Proposition 14.41 and 
Theorem 1121 let = Q^{^r) and set 

E{x,u) = ^e,o,x{u), defined on K x ^ M". 

To establish assumption (i) of Theorem 14.51 it suffices to note that by virtue of condition (iii) in 
Theorem 14.21 one has that for x G B^{xo,r), 

B,{xo,2r) C B,{x,3r) C E{x,Q,). 

Assumption (ii) in Theorem 14.51 is a direct consequence of condition (ii) in Theorem 14.21 with 
ai = 16. By the classical connectivity result of Chow we see that E{x, u) satisfies assumption 
(iii), with a function 7, piecewise expressed as exponential mappings of vector fields of e— degree 
one. Let us denote the required vector fields. With this choice of path, it is known (see 

for example [ini Lemma 2.2] or [9l pp 99-101]) that x — )• ^{x,u,t) is a path, with Jacobian 
determinant 

1 + 'il^{x,u,t), 



(ieiTT— (x, u, t) 
ox 



for a suitable function ^p{x, u, t) satisfying 

\ip{x,u,t)\ < cr, on K X X [0,cr]. 

Moreover the constant c only depends on the Lipschitz constant of the vector fields {Xf)i^i^. Hence 
in our setting it can be chosen independently of e. Consequently also condition (iv) is satisfied. 

□ 
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